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Abstract We study the kinetics of confined carrier-phonon system in a quantum 
dot under fast optical driving and discuss the resulting limitations to 
fast coherent control over the quantum state in such systems. 



1. Introduction 

Unlike natural atoms, semiconductor quantum dots (QDs) always 
form part of a macroscopic crystal. The interaction with the quasi- 
continuum of lattice degrees of freedom (phonons) constitutes an inher- 
ent feature of these nanometer-size systems and cannot be neglected in 
any realistic modeling of QD properties, especially when the coherence of 
confined carriers is of importance. The understanding of the decisive role 
played by the QD coupling to lattice modes has increased recently due to 
both experimental and theoretical study (spectrum reconstruction [1—4], 
relaxation [5-10], phonon replicas and phonon-assisted transitions [11- 
16], phonon-induced pure dephasing upon ultrafast excitation [17-22]). 
The phonon-induced decoherence seems to be crucial for any quantum 
information processing application and for any nanotechnological device 
relying on quantum coherence of confined carriers [23, 24]. 

There are three major mechanisms of carrier-phonon interaction [25]: 
(1) Coulomb interaction with the lattice polarization induced by the rela- 
tive shift of the positive and negative sub-lattices of the polar compound, 
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described upon quantization by longitudinal optical (LO) phonons; (2) 
deformation potential coupling describing the band shifts due to lattice 
compression, i.e. longitudinal acoustical (LA) phonons; (3) Coulomb in- 
teraction with piezoelectric field generated by shear crystal deformation 
(transversal acoustical, TA, phonons). The lattermost effect is weak in 
InAs/GaAs systems but may be of more importance for the properties 
e.g. of GaN dots [26-28]. 

Due to the carrier-phonon interaction, any change in the carrier sub- 
system must be accompanied by the appropriate modification of the 
lattice state. An example of this effect may be creation of the exciton- 
polaron states: excitons accompanied by the lattice polarization (coher- 
ent LO phonon field) with energy shifted down by a few meV [3, 22] . 
The interaction with acoustical branches leads to creation of a similar 
deformation field which, although much less pronounced in terms of en- 
ergy shift, strongly influence the dynamical properties of the interacting 
system due to gapless spectrum and low characteristic frequencies of the 
acoustical phonons. 

The lattice response to any manipulation performed on the charge 
distribution confined in a QD leads to discrepancy between the desired 
state of the quantum confined system and the actual one. This effect 
may be stronger than any decoherence process in an undriven system. 
In terms of the quantum information processing schemes, such a dis- 
crepancy is manifested by the loss of fidelity of the quantum operation. 
In quantum nano-technology applications, it will reduce the efficiency of 
nano-devices. 

In the present paper we discuss the carrier-lattice kinetics induced 
by optical excitation of a confined exciton. The paper is organized as 
follows: In the next section we present the formal model of the exciton 
confined in a QD. The section 1.3 describes the system kinetics in re- 
sponse to an abrupt change of the charge distribution. In the section 
1.4 we discuss the trade-off between the dynamically induced error and 
other processes limiting the coherence time of a quantum state. The 
final section contains concluding discussion. 

2. The model 

We consider the Hamiltonian describing a single exciton interacting 
with phonons, 
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Electron mass 


m c 


0.067m 


Hole mass 


m h 


0.38m 


Static dielectric constant 


£s 


13.2 


Effective dielectric constant 


i 


62.6 


Optical phonon energy 




36 meV 


Longitudinal sound speed 


c 


5150 m/s 


Deformation potential for electrons 


(To 


6.7 eV 


Deformation potential for holes 




-2.7 eV 


Unit crystal cell volume 


V 


0.044 nm 3 


LO phonon dispersion parameter 


P 


0.03 meV-nm 2 



Table 1. The material parameters used in the calculations (partly after Refs. [29, 
30]). 



where 6 S) k, b\ k are bosonic annihilation and creation operators for LO 
(s = o) and LA (s = a) phonons with quasi-momentum k. The cor- 
responding frequencies within the effectively coupled wavevector range 
may be modeled by w 0j k = — Q — fik 2 and w a> k = ck. The coupling 
constants for the LO and LA phonon branches are given by 



and 

n(a) 



(i) 

$;(r c ,r h ) [aj^ - a h e ikrh ) $ n /(r c , r h )d 3 r e d 3 r h , 



y 2gvc 

where r e , denote the coordinates of the electron and hole, respectively, 
and $ n (r e ,rh) is the exciton wavefunction. The other elements of the 
notation are described in Table 1, along with values (corresponding to 
InAs/GaAs system) used in the calculations. 

Numerical diagonalization of the interacting electron-hole system in 
parabolic confinement, under assumption that non-interacting electron 
and hole would have the same wavefunctions, leads to the spectrum 
shown in Fig. la (M is the conserved total angular momentum). The 
dominant contribution to the lowest excited states of the exciton comes 
from the excited hole states, while the electron wavefunction is only 
slightly modified. The corresponding electron and hole distributions are 
shown in Fig. lb-d. 

The exciton-phonon coupling functions may be calculated using the 
wavefunctions found numerically. Fig. 2 (left panels) shows the results 
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Figure 1. The spectrum of the exciton in QD (upper left) and the electron (solid 
lines) and hole (dashed lines) probability densities for the three lowest exciton states, 
compared to the ground state probability density for a noninteracting particle (dotted 
lines). Here the noninteracting particle wavefunction width is l± = 4.9 nm in-plane 
and l z = 2 nm in the growth direction. 

for coupling between the ground state and a few lowest states, averaged 
over angles. The coupling to LO phonons is much stronger than to LA 
phonons. In the case of LO phonons, the coupling strength increases for 
excited states due to less charge cancellation between the electron and 
hole in these states (cf. Fig. lb-d). 

3. System kinetics after an ultrafast pulse 

Let us study the kinetics of the system after an excitation performed 
by an extremely short (formally infinitely short) laser pulse. The in- 
formation on the dynamics of the interacting carrier-lattice system is 
contained in the exciton single-particle causal Green function, 

Gn,n'(*) = -i(T{a n (t)al(0)}). 

The averaging (...) is the temperature-dependent averaging with respect 
to phonon degrees of freedom and the vacuum of exciton (cf. Ref. [31]; 
it corresponds to the case when the grand canonical averaging sector 
without exciton, vacuum, is energetically distant, ~ 1 eV, from the next 
sectors). For t > this Green function coincides (up to a constant 
factor) with the correlation function (a n (t)a n ,(0)} which, for n\ = ri2 
may be interpreted as a measure of integrity of the excitonic state: it 
corresponds to the overlap of the state at time t with this state at the 
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Figure 2. Coupling constants for LA (a) and LO (b) phonon modes between the 
ground state and the three lowest states and the resulting structure of the imaginary 
part of the mass operator (c,d). 



initial moment t = 0. The Fourier transform of the correlation function, 
I n)n i{uj) = J^° oo (a n (t)al l ,(0))e lult dt, is usually called the spectral density 
[32, 33], and it can be expressed by the imaginary part of the causal 
Green function: 

ImGVH = ~\ (l + e-/^ T )~ 1 /„ in ,H, (2) 

where G n , n ,(u) = JZo G n , n/ {t)e iu,t dt. 

The equation of motion for the causal Green function may be rewritten 
as a Dyson-type equation, 

(u - e n )G n , n ,{u)) - M n ^t{u)G n ,i jn ,(u) = (3) 
with the mass operator 

Mn,A») = ^ E ^(k)*£U-k) 

k,s,n",n"' 

x J dw'G n >/ ;n >//(uj + lu')D°(\l,s,uj'), 

where we have restricted ourselves to single-phonon processes by replac- 
ing the vertex function by its zeroth order approximation and using the 
free phonon Green function (LO phonon damping may be easily included 

later: see below). Since F^,(k) ~ g s , then |M ni „/(u;)| 2 ~ g 2 s . Hence, up 
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to gj one has from Eq. (3), 
and 

M nn {u) = A n (uj) - ij(uj) 



uj- e n - M n)Tl (uj) 



(4) 
(5) 



2vriV 



J2\ F nn>(V\ 2 [ ^G n ^u, W)£>^( S ,k,u/). 

n',k J 



For the real and the imaginary part of the mass operator we obtain 
form (5) 



iV 



(6) 



k,s,n, 



+ 



(1 + n Sik )(a> - e w / - A n /(o; - ^ S;k ) - u; Sik ) 
[w - e n > - A n /(w - w S)k ) - w Sik ] 2 + 7^,(0; - w s (k)) 

ra S;k (o; - e n / - A n /(o; + a^k) + ^,k) 
[u> - e n > - A n / (a; + w Sjk ) + w Sjk ] 2 + 7 2 , (m + w s>k ) 



and 

7nH 



^ E i*$ooi 



E 

k,s,n' 



(7) 



(1 + n S;k )7 re /(u; - u^k) 



[w - e n / - A„/(w - w Sjk ) - w S)k ] 2 + 7 2 ,(w - w s (k)) 



A n /(w + ^ s , k ) +w Sjk ] 2 + 7 2 , (a; + w Sjk ) 



In this equation the first term determines the energy shift due to dress- 
ing with LO phonons, while the second one corresponds to LA phonons. 
The former dominates energetically: the energy shift is mainly due to 
the interaction of exciton with optical phonons (dressing of exciton with 
LO phonons, creating together an exciton-polaron) . The latter is small 
and its energetical effect may be safely neglected but it contributes con- 
siderably to the system kinetics. 

To the lowest order, one may neglect A n / and 7„/ in the right-hand- 
side of (6) and, neglecting the acoustic phonon term, write the equation 
for the exciton-polaron energy levels E n as the poles of the Green func- 
tion (4), E n -e n - A n (E n ) = 0, i.e. 



1 

N 



1 + ™ ,k 



k,n 



En 



+ 



n o ,k 



$7 E n — e n + £1 



0. (8) 
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At kv,T <C the above equation is equivalent to that found by Davydov 
diagonalization of the Frohlich Hamiltonian for exciton [34]. The Eq. (8) 
for n = gives the ground state energy shift Ao ~ —5 meV for the QD 
with parameters listed in Table 1 (see Ref. [22]). 

The imaginary part of the mass operator is given to the leading order 
by the equation (putting 7 n / = in the rhs. of Eq. (7)): 



7nM = (9) 

Jf E {l^SflOl 2 K 1 + n oM)$(u ~ E n , - fi k ) + n oM 5(u; - E n , + fl k )] 

k,n 

+\F® (k)| 2 [(1 + n a , k )5( W - £7 n , - cfc) + n a , k 5(w - £ n , + cfc)]} , 



where we use the fact that the equation u— e n '— A n /(o;±a; S)k )±a; Sik = is 
solved by u = E n > fw s ^ and neglected higher-order corrections resulting 
from resolving the Dirac 5. The first term in Eq. (9) describes the energy 
transfer to the LO phonon sea, while the second one corresponds to the 
energy transfer form gradually dressing exciton to the LA phonon sea. 

The form of 70 (uj) for our model, including a few lowest states, ob- 
tained using the numerical wavefunctions is shown for various tempera- 
tures in Fig 2c, d. 

Using (4), one finds 



imo n n w — — j— 



(10) 



In the following, we will focus on the ground exciton state, n = 0. Usu- 
ally, this state is long-living at low temperatures: the broadening of the 
corresponding spectral line, related to radiative lifetime and thermally 
activated processes, does not exceed 0.1 meV for T < 100 K [35, 17]. 

Let us note that the term Jq(u>) in the denominator of (10) is im- 
portant only near Eq, where the other term vanishes (otherwise, it is a 
higher-order correction). Therefore, it may be approximated by a con- 
stant 70 = 7o(-Eb)- in a similar manner, the correction Aq(uj) may be 
replaced by Aq(uj) Ao(-Eb) + — Eo)(dAo(uj)/duj) u =E an d combined 
with eo to give Eq, in accordance with (8). 

Thus, one may write for Iq(lo) = — (l/7r)ImC?oo(^) 

I (u) = (11) 

z -l 1 70/2 



vr (oj - EqY + 7o 2 /4 
-r- . ( s ). s |2 ( n s,k + - E n - w s k ) + n S;k S(uj - E n + w S)k ) 

+ „?' F °» (k)l 
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Figure 3. Left: Spectral density vs energy (inset shows the shape of the LO phonon 
replica) Right: The corresponding evolution of the correlation function for two tem- 
peratures, as shown (Insets show the evolution on longer timescales). The calculation 
is limited to the two lowest exciton states. 



where 

dAjuj) 1 
Z - 1 ~-fo-\"=*>~ 1 +N 2^ 

k,s,n' 

The Fig. 3a,c shows the result. 

At T = 0, 70 (w) = at lo = Eq (neglecting recombination process) 
and this point is a well defined pole of the causal Green function. It 
corresponds to a quasiparticle: the dressed exciton (i.e. the exciton- 
polaron, if neglecting LA phonons). For T > 0, this spectral density has 
a Lorentzian shape around lo ~ E , with width 70 . Our model is unable 
to quantitatively account for the broadening found in experiment. In 
fact, to our knowledge, its origin has not been explained so far; we will 
use the experimental values of 70 corresponding to 630 ps and 170 ps 
exciton lifetime (including radiative porocesses and thermally induced 
transitions) at T = K and T = 25 K, respectively [17]. 

Apart from the central peak, the spectral density shows acoustic and 
optical phonon sidebands. To the former, only the ground state itself 
contributes: the magnitude of the features visible in 7(0;) around excited 
states (Fig. 2c, d) is much smaller than the energy distance between these 
states. Therefore, as far as the interaction with acoustical phonons is 
concerned, a single level (independent boson) model supplemented by 
central line broadening is very accurate. 



[1 + 2n s , k ] 



Fast control in quantum dots 



9 



The features resulting from the LO phonon branch behave in a dif- 
ferent manner: the contribution from nonadiabatic coupling to excited 
states is considerable. This results from the fact that the coupling to 
LO phonons is much stronger in general, the energies of LO phonons are 
comparable to the exciton energy spacing and the coupling to excited 
levels is stronger than to the ground state itself due to partial charge 
cancellation in the latter. In this case, the contribution from the higher 
states dominates over the ground state contribution and the independent 
boson approach is not valid. 

The time evolution of the dressing process is described by the inverse 
Fourier transform of the correlation function given by (2). Taking into 
account that for u in the energy sector of exciton (of order of eV) u » 
&bT, one has 

1 f°° 

In{t) = — MmG n , n ( W )e-^. 
vr J-oc 

The time-dependent correlation function may be obtained from (11), 
Io(t) = 

z -l e ~i(E -i~f /2)t 

(n k + l) e - i ( s "+^,k)i n ^ e -i{E n -oj sM )t 
{AE n0 + co s , k ) 2 + 7o 2 + (AE n0 - oj s ^Y + 7 g J ' 

where AE^o = E n — E$. The result, calculated numerically for various 
temperatures, is plotted in Fig. 3b, d. 

The smooth sidebands of I(cu) correspond to the initial correlation 
decay on 1 ps timescale. The LO phonon peaks lead to slowly decaying 
fast oscillations manifesting phonon beats with frequencies correspond- 
ing to the energy differences E n — Eq — Q. Without including anhar- 
monic effects, these oscillations would decay very slowly due to the weak 
LO phonon dispersion. When the anharmonic phonon damping (corre- 
sponding to tlo = 9 ps decay of an LO phonon [36, 37]) is included 
by substituting w ,k - ► <^o,k ± &7lo> the LO phonon beats decay with 
characteristic time slightly below tlo, which results from the joint effect 
of damping and dispersion. The long-time dynamics is governed by the 
Lorentzian feature around lo = Eq and shows an exponential decay (as 
assumed beforehand). 

4. System response for finite-duration pulses 

In the previous section we have shown that any change in the carrier 
distribution in a QD is followed by lattice relaxation (dressing) processes 
which lead to some loss of the quantum coherence. To preserve the coher- 
ence while operating on the carrier states, the action must be adiabatic 
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with respect to lattice timescales, i.e. phonon periods. Slowing down 
the dynamics leads, however, to increasing effect of other decoherence 
mechanisms, like radiative recombination or thermally activated transi- 
tions to higher states, as the decoherence caused by such effects increases 
linearly with time. 

In order to study the interplay of these decoherence effects we have 
studied the evolution of an exciton coupled to acoustic phonons under 
optical excitation by a finite-duration laser pulse [23]. The effect of 
coupling to phonons was quantified in terms of the error of a coherent 
operation, 8 = 1 — F, where the fidelity is defined as an overlap between 
the actual state and the desired one (the latter was taken to be the 
dressed final state, obtained by adiabatic operation), F = (cp\p\(j)}, where 
\4>) is the ideal final state and p is the final density matrix. 

The density matrix for the final state is calculated using second order 
expansion for the evolution operator of the total system and tracing over 
the lattice degrees of freedom (see [23, 38] for details). For timescales 
relevant here, it is sufficient to consider acoustic phonons and only one 
(ground) exciton state. The error averaged over initial states can be 
represented as the overlap of two functions 

6 = J ^R{u)S{u), 
where R(u>) is the spectral density of the reservoir 

R(u) = £ [Sfa - w)(n k + 1) + 5(w k + co)n k ] |F(k)| 2 . 

k 

The function S(cu) represents spectral characteristics of the system and 
is given by 

5H = i(|F_H| 2 + |F + M| 2 ), |F ± H| 2 «aV |r| (^*) j (12) 

for a Gaussian driving pulse f(t) = a/(\/27rTg)e~ 2 (*/ r s) ; rotating the 
state by the angle a on the Bloch sphere. 

For the coupling (1) to LA phonons, one finds at T = for uj <C c/J 1 

RM^BtfS, Ro= ( -^f, (13) 

and we obtain the averaged error 5 = a 2 RoTg 2 /3 (taking into account 
the upper cut-off, the error will be finite even for an infinitely fast gate 
(see Fig. 4). 
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Hence, this non-Markovian error increases as the speed of gate in- 
creases. This could result in obtaining arbitrarily low error by choosing 
suitably low speed of gates. However, if our system is also subjected to 
other type of noise this becomes impossible. Indeed, if we assume that 
our system undergoes additional amplitude damping with rate 7m, the 
total error per gate is 

= — 2~ + 7MT g , 7nM = -a R , 7m = — , (14) 

where r r is the characteristic time of exponential decay (Markovian er- 
ror). As a result, the overall error is unavoidable and optimization is 
needed. The above formulas lead to 

<W = ^7nM) 1/3 = I (^) 1/3 = « 2/3 0.0035, 

for 

r g = (2^) V3 = Q« 2 «or r ) 1/3 = a 2 / 3 1.5ps, 
where we have used GaAs material parameters (Table 1). 
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Figure 4- Combined Markovian and non-Markovian error for a a = n/2 rotation on 
a qubit implemented as a confined exciton in a InAs/GaAs quantum dot, for T = 
(solid lines) and T = 10 K (dashed lines), for two dot sizes (dot height is 20% of its 
diameter). The Markovian decoherence times are inferred from the experimental data 
[17]. Inset: Spectral density of the phonon reservoir R(w) at these two temperatures 
and the gate profile S(ui) for a = n/2. 
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The exact solution within the proposed model, taking into account the 
cut-off and anisotropy of the dot shape and allowing finite temperatures, 
is shown in Fig. 4. The size-dependent cut-off is reflected by a shift of 
the optimal parameters for the two dot sizes: larger dots admit faster 
gates and lead to lower error. Interestingly, the trade-off becomes more 
apparent at nonzero temperature. 

It should be noted that these optimal times are longer than the limits 
imposed by level separation [39-43]. Thus, the non-Markovian reservoir 
effects (dressing) seem to be the essential limitation to the gate speed. 

5. Conclusions 

We have shown that an action performed on the carrier system (e.g. 
exciton) confined in a quantum dot must be accompanied by an appro- 
priate reconstruction of the lattice state. Such a dressing process takes 
place spontaneously if the charge distribution is changed on times short 
compared to timescales of the lattice dynamics (phonon oscillation peri- 
ods). This spontaneous process destroys the quantum coherence of the 
system and may preclude quantum information processing or coherent 
nanotechnology applications. 

On the other hand, sufficiently adiabatic operation would take long 
time compared to the exciton lifetime and other exponential decay pro- 
cesses, again leading to large coherence loss. Therefore, optimal operat- 
ing conditions for maximum preservation of coherence must be searched 
for. 

We have evaluated the minimal error for typical parameters, showing 
that the optimal gating time (~ 1 ps) lies within the current experimen- 
tal possibility. This optimal time sets up the limit beyond which further 
gate speed-up is unfavorable. It is remarkable that minimal error 5 we 
found (0.004 — 0.01), although not extremely high, is still of 1-2 orders 
of magnitude higher than the one admitted by fault-tolerant schemes 
known so far (~ 10~ 5 ). However, possible improvements of the latter 
schemes cannot be excluded. 
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